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We report on theoretical investigations of scattering asymmetry vs. incidence of carriers through 
exchange barriers and magnetic tunnel junctions made of semiconductors involving spin-orbit inter¬ 
action. By an analytical 2x2 spin model, we show that, when Dresselhaus interaction is included in 
the conduction band of antiparallel magnetized electrodes, the electrons can undergo a large differ¬ 
ence of transmission depending on the sign of their incident in-plane wavevector. In particular, the 
transmission is fully quenched at some points of the Brillouin zone for specific in-plane wavevectors 
and not for the opposite. Moreover, it is universally scaled by a unique function independent of the 
spin-orbit strength. This particular feature is reproduced by a 14 x 14 band k • p model showing, 
in addition, corresponding effects in the valence band and highlighting the robustness of the effect, 
which even persists for a single magnetic electrode. Upon tunneling, electrons undergo an asym¬ 
metrical deflection which results in the occurrence of a transverse current, giving rise to a so-called 
Tunnel Hall Effect. 

PACS numbers: 72.25.Mk, 75.70.Tj, 75.76.-|-j 

The interplay between particle spin and orbital mo¬ 
tion is currently at the basis of new functionalities. 

These generally require efficient spin-current injection at 
magnetic-nonmagnetic interfaces, efficient spin-transfer 
torque, and possibly efficient spin Hall effect with heavy 
materials mi for magnetic commutation. Spin-orbit in¬ 
teraction (SOI) at an interface with broken inversion 
symmetry may lead to the observation of Rashba-split 
states [SHE] which can be used to convert a perpendicu¬ 
lar spin-current into a lateral charge current by inverse 
Edelstein effect {I EE) 0 HD]. In that context, inves¬ 
tigations of SOI in solids and at interfaces is of prime 
importance for basic physics and today’s technology. 

In this paper, we study unconventional quantum ef¬ 
fects resulting in giant transport asymmetry for electrons 
or holes in structures composed of exchange and spin- 
orbit-split electrodes made of HI-V semiconductors with 
antiparallel (AP) magnetizations, possibly separated by 
thin tunnel barriers. The symmetry of the structure al¬ 
lows a transmission difference vs. carrier incidence near 
interfaces with respect to the reflection plane defined by 
the magnetization and the surface normal. This quantum 
process departs from the effect of a beam deviation by the 
Lorentz force due to the action of a local magnetic field in 
the barrier m and from spin-filtering effects occurring 
in non centro-symmetric structures mis]. Unlike the 
latter, the effect we propose requires the simultaneous 
action of both in-plane and out-of plane spin-orbit fields 
for promoting transport spin asymmetry. In order to ad¬ 
dress the issue in a simple way, we first consider a hetero¬ 
junction made of two identical magnetic semiconductors 
of zinc-blende symmetry, with opposite in-plane mag¬ 


netizations: this structure (Fig. 1) constitutes an ideal 
exchange step and a paradigm for exchange-engineered 
heterostructures, similarly to symmetrical spin valves in 
giant magnetoresistance mug. Indeed, due to the ax¬ 
ial character of the magnetization, the AP configura¬ 
tion breaks the symmetry with respect to the reflection 
plane (in Fig. 1 the reflection plane is the xz plane), 
and also some possible rotation and time conjugation in¬ 
variances existing in the parallel (PA) magnetic arrange¬ 
ment m- The result is that two waves with opposite 
in-plane wavevectors ky may be differently transmitted 
in amplitude. 

We first consider Dresselhaus interaction in the con¬ 
duction band of bulk materials ini- Hereafter, we re¬ 
fer the structure to the x, y, z cubic axes (unit vectors 
X, y, z) and assume that electron transport occurs along 
the z axis, whereas the magnetization lies along x. We 
study the transmission asymmetry when the wavevector 
component along y is changed from ^ to —Electrons 
are injected from the first conduction band of material I 
to the left (e = 1) into the first conduction band of ma¬ 
terial II to the right (e = —1). Then, the relevant 2x2 
Hamiltonians respectively write: 

Hiji = 7c + i^)l + wm ■ a + (yy) • a 

^ / 7c(fc^-7^^* -iU^k'^ + ew V , 

y ij^k^ + ew 7c(fc^ + + 7^^^ / 

where (0, k) is the electron wavevector; I is the iden¬ 
tity matrix, 7 c accounts for the conduction effective mass, 
m is the unit magnetization vector, 2w the exchange 
splitting (assumed to be positive), a the Pauli opera¬ 
tor, and X = the D’yakonov-Perel’ (DP) 
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Figure 1 . Scheme of transmission process at an exchange-SO/ 
step with AP magnetizations M and —M along x. The prop¬ 
agation direction of carriers (straight arrow) is along 2 with 
propagative wavevector fci whereas the in-plane incident com¬ 
ponent (heavy line) or —^ (dashed line) is along y, xyz 
forms a direct frame. The dash-dot curve denotes the evanes¬ 
cent waves, either reflected or transmitted. Carriers with 
in-plane wavevector component are more easily transmitted 
than those carrying —^.(Top right inset): Energy profile of 
the exchange step; E is the longitudinal kinetic energy along 
0 and 2 w Is the exchange splitting in the magnetic materials. 


internal field responsible for the spin splitting mM- 
For the subsequent discussion, we introduce the tensor 
7 = i'jiSij) which characterizes the DP-field strength, 
with 7a; = 7y = 7, 72 = 7; and Sij the Kronecker symbol. 
We consider the two cases 7 = 7 and 7 = 0 , switching on 
and off the diagonal perturbation. 

The two energies in the exchange and spin-orbit-split 
subbands are given by £i = 7^ [kf + — w and £2 = 

7 c(fci + ^^) + w, where ki (^2) is the z-component of 
the wavevector in the lower (upper) subband. These 
expressions are correct up to first order in 7 provided 
|7^^fc/u;| << 1 and |7^A:^/iu| << 1, where k = ki or 
k = k2- The respective eigenvectors write: 

W.i (^, ki) = [1 - 2ei^lkl,-e (1 - 2Ji^ki)] /v^, (2) 

Ue ,2 (C, ^ 2 ) = [1 - 2€ifikl, e (1 -f 2^^fc2)] /V2, (3) 

where /i = ^^/{ 2 w) and Ji = j^/{ 2 w). Note that the 
norm of (^ = 1 or 2) only involves even powers of 
^ likewise the direct overlap |(ue_^|u_£7)|^ between in¬ 
coming and outgoing states, so that no transmission 
asymmetry can be expected in usual tunneling models, 
e.g. based on interface density of states [MSI]- The 
asymmetry appears in full-quantum treatments involving 
matching conditions at interfaces and may be correctly 
described by embedding methods |22| . 

The corresponding wavefunctions in Regions I and II 
can be written in a compact form: 

^'7 (z) = aui.2 (^, ^2) -k / 3 uyi (^, ki) 


+Aui, 2 (e, -k2) + Rui,i (e, -fci) 


-ikiz 
^ 7 


^7/W =Cu_i.i(^,A:i) e*'=i^-ki5u_i,2(C,/c2) 

( 4 ) 

where the a and fj (resp. A and B) amplitudes stand 
for incident waves (resp. reflected waves) in Region I, 
and C and D for transmitted waves in Region II. Be¬ 
cause k\\ is conserved in the transport process, we are 
dealing with states with the same longitudinal kinetic en¬ 
ergy E along z and a total kinetic energy £ = E + 7 c^^. 
The proper matching conditions are the continuity of 
the wavefunction and of the current wave J'I'777 = 
(l/?i) {dHjji/dk)'^jjj because Hjjj contains no more 
than quadratic k terms | 23 tl? 7 ] and because 7/7C is con¬ 
tinuous. The average transmission coefficient T (^, ^1,^2) 
upon positive and negative incidences is related to the 
amplitude of the transmitted wave C (^, ^1,^2) calculated 
with the initial conditions a = 0 and j 3 = \ through: 


T{^,ki,k 2 ) 


|C(e,fci,A:2)|V|C'(-e,fci,A:2)|' ,,, 


and we define the transmission asymmetry as: 


A{^,ki,k 2 ) 


|C(g,fcl,fc2)|^-|C(-C;fcl,fc2)|" 

|C(^,fci,fc2)|V|C(-C,A:i,fc2)r 


( 6 ) 

It can be checked that, when 7 = 0 , (^, fci, ^2) van¬ 

ishes if a and /3 are real, which is a non trivial result. 
The transmission of a pure up-spin incident electron into 
a pure down-spin state is only possible under oblique in¬ 
cidence via SOI which introduces off-diagonal matrix el¬ 
ements. Moreover, a non-vanishing diagonal part of SOI 
is necessary to obtain a non-zero asymmetry although 
the z component of the DP field along z does not depend 
on the sign of k\\ | 28 | . Then, from now on, we take 7 = 7. 
The wavevector ki in the lower subband has to be real 
so that we can define K = ki > 0 . We introduce the pa¬ 
rameter A > 0 with ^2 = iXK, the reduced longitudinal 
energy g = E/w = (l — A^) / (l-|-A^), and the incidence 
parameter t = ^/K. One obtains: 


[tf/A-)(3V-l)+2A(V-l)] _ 

\ W J (A-f) 

( 7 ) 

From Eq. it can be checked that A (^, fci, ^2) = 0 if 
A is purely imaginary; the asymmetry appears when the 
lower-energy band carries a propagative state whereas 
the upper one acts as a barrier sustaining an evanescent 
state. Transport is then described in a two-/c-channel 
model, a propagative channel (fci) and an evanescent 
channel {k2). One obtains: 

T {t, g) = (1 -k gf [Ig^ {1 - g) + + g){2g - 1)^] , 

Tc 

( 8 ) 
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Figure 2. (a) Universal asymmetry coefficient A vs. reduced 
energy r] = E/w obtained for different values of t = ^/K 
[t = 0.01 (black; circles), t = 0.5 (blue; squares), t = 1 (red; 
stars), and t = 2 (purple; triangles) by 2 -band analytical (full 
line) and numerical (symbols) calculations. Two-dimensional 
map of the transmission coefficient T in 2 x 2 (b) and 14 x 14 
(c) k ■ p band models for the exchange-SO/ step schematized 
in Fig. 1; the parameters are: exchange energy 2w = 0.3 eV, 
total kinetic energy £ = 0.08 eV counted from the middle of 
the conduction step, and 7 = — 24 eV A®; band parameters of 
the 14-band k • p model taken from Ref. |30| 


and 


A{t,r]) 


4 ^ 771/1 — 77 ^ (2r] — 1 ) 

4772 (1 — 77 ) -I- (1 -I- 77 ) (277 — 1)^ 


( 9 ) 


Equation emphasizes the increase of T (t, 77 ) with t 
and 7 . The range of validity defined above can be written 
\ « 1 , a condition easily fulfilled. 

The analytical asymmetry A is plotted in Fig. 2a for 
several values of t (full lines), where the symbols refer 
to the 2 x 2 numerical calculations, showing an excellent 
agreement. It can be seen that the curves related to t 
and t = 1 /t are located at almost symmetrical positions 
with respect to the t = 1 curve. They admit four zeros 
in the energy range considered: i) two at the two ends 
of the energy step when either the propagative or the 


evanescent state disappears and ii) one in the middle of 
the energy barrier and one for r] = 3 ji which is particular 
to Dresselhaus interaction. It is a remarkable result that 
A {t, rj) does not depend either on the material parame¬ 
ters or on the sign of 7, thus conferring to a universal 
character. Reversing the magnetization (changing w into 
—w) makes transport occur in the k2 channel leading to 
a change of A{t,rj) into —A{t,r]) [ 29 ]. An important 
result is that A is positive when (m, k) forms a direct 
frame and negative otherwise. Another striking feature is 
that an arbitrarily small perturbation is able to produce 
a 100% transport asymmetry i.e., a total quenching of 
transmission. Fig. 2 b- 2 c display the 2 -dimensional map 
of the electron transmission at a given total energy in 
the reciprocal space calculated using both a 2 x 2 ef¬ 
fective Hamiltonian (Fig. 2 b) and a full 14 x 14 band 
k-p treatment (Fig. 2 c) involving odd-potential coupling 
terms P and A | 30 llH 2 ] . These calculations are based 
on the multiband transfer matrix technique [231 EH Eg. 
We have checked that transport asymmetry also arises 
for a tunnel junction where a thin tunneling barrier is in¬ 
serted between the two magnetic layers. Tailoring more 
complicated structures which involve resonant tunneling 
effects allows one to obtain much higher transmission (a 
fraction of unity) while keeping the same magnitude for 

A Eg. 

We can calculate the transmitted current, J [t, 77 ] = 
Jj [iF//(z)] -I- J_.c [iF//(z)], originating from incident 
waves of equal amplitude with opposite fey. To the lowest 
order in 7, we find 

V] = -(1 -b 77)^/^ T (t, 77) [A it, ri) ty+z] 

( 10 ) 

Thus, the asymmetrical transmission gives rise to a trans¬ 
verse momentum and then to a tunneling surface current 
(per unit length) jy = Jy x imfp {^mfp is the electron 
mean free path) which can lead to an anomalous Tunnel 
Hall Effect (THE) under steady state regime. This effect 
could be experimentally investigated at a scale where the 
thickness of the channel collecting the current is compa¬ 
rable to imfp, *-e-> not exceeding a few nm | 35 | . The ra¬ 
tio of the (surface) transverse to the longitudinal current 
jy \j,v]/Jz [i-i'n] = tA{t,r]) imfp then defines the THE 
length in the spirit of a recent work dealing with lEE 
phenomenon HED]. To gain some numerical insight, an 
incident beam in Region I with a given angular disper¬ 
sion with respect to the z axis gives rise, after angular 
averaging at a fixed given energy |3g, to a THE length 
tAim fv of the order of im fn for a beam deviation as large 
as 45 “ Eg. 

What happens for the valence bands simply described 
within the 6-band Luttinger k • p effective Hamilto¬ 
nian EH? The results are shown in Fig. 3 ; we have 
checked that the 14 -band model provides similar data 
with P = 0 and A = 0 , surprisingly showing that the 
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absence of inversion symmetry is not a key feature in 
the valence band. The lower curve displays the asym¬ 
metry A vs. hole energy £ in the case of a 3-nm-thick 
tunnel barrier. The energy range covers the valence spin 
subbands, namely, starting from the highest energy, the 
up-spin heavy (light)-hole band HH ^ {LH f), the down- 
spin light (heavy)-hole band LH (HH |), and the up 
(down)-spin split-off band SO t(5'0 1). We refer to 
points (1) to (6) marked by vertical arrows in the fol¬ 
lowing discussion. Here, the energy of the HH f [HH J,] 
maximum, corresponds to 0.15 eV [—0.15 eV], the en¬ 
ergy origin being taken at the top of the valence band 
of the non-magnetic material, and is indicated by point 
(1) [(4)]. Correspondingly, one observes an almost fully 
negative transmission asymmetry in this energy range for 
predominant majority up-spin injection as far as HH j, 
does not contribute to the current. At more negative en¬ 
ergy \£ < — 0.15 eV: point (4)], a sign change of A occurs 
at the onset of HH j, (in the upper left inset, see the 
step in the transmission coefficient, which reaches almost 
+50%); A remains positive after crossing SO f [point (5)[ 
before turning negative again once crossing 5'0 j, [point 
(6)[. Note that A changes sign two times at characteris¬ 
tic energy points corresponding to a sign change of the 
injected particle spin. We have performed similar calcu¬ 
lation for a simple contact (d = 0; right upper inset in 
Fig. 3, black curve). Remarkably, A, although smaller, 
keeps the same trends as for the 3-nm tunnel junction, 
except for a change of sign, showing a subtle dependence 
of the exchange coupling on the barrier thickness. With¬ 
out tunnel junctions, A abruptly disappears as soon as 
SO j, contributes to tunneling [circle region[ i.e., when 
evanescent states disappear. In the case of tunnel junc¬ 
tion, A, although small, subsists in this energy range and 
this should be related to the evanescent character of the 
wavefunction in the barrier. 

To gain physical insight in the case of the valence 
bands, let us start from a semiconductor toy model with 
conduction and valence states described by Kane’s S, 
X, Y, and Z real cubic harmonics (no SOI) [3H]- As¬ 
sume that the exchange-split semiconductors are sep¬ 
arated by a barrier involving SOI, hereafter intro¬ 
duced as a perturbation. Then, the relevant wave- 
functions in the barrier are respectively associated to 
the wavevectors ± iXKz = K {tyzL iXz) to the left 
(+) and right (—). To second order, k • p perturba¬ 
tion shows that only the Y and Z valence wavefunc- 
tions become coupled through the conduction state, lead¬ 
ing to two energy subsets with eigenstates +*AF + tZ 
(no energy change) and tY ± iXZ with an energy shift 
in the gap equal to (A^ — {KP)^ /Eq, with = 

| 27 c/?i (S'jpzjZ) p. In the AP state, (—iXY + tZ) f to 
the left [resp. {tY + iXZ) t[ and {iXY + tZ) j, to the 
right [resp. {tY — iXZ) ![ belong to the same energy 
shell, resulting in efficient matching at particular inci¬ 
dence t. These eigenstates are carrying orbital momen- 
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Figure 3. {Bottom)-. Transmission asymmetry A vs. total 
energy £ for a magnetic tunnel junction in the AP state. 
The parameters are: 2w — 0.3 eV, parallel wavevector ^ = 
0.2 nm“^, barrier thickness d = 3 nm, and barrier height 0.6 
eV. The energy zero corresponds to the non-magnetic upper- 
valence-band maximum. {Upper left inset)-, transmission T 
calculated in the AP state; {Upper right inset)-. Asymmetry 
A without tunnel barrier (d = 0 ; black) compared to the case 
of the tunnel junction (dotted red). 


turn along x, namely (Ax) proportional to — At, whereas 
{Ly) = {Lz) = 0. Branching L • S SOI into the barrier 
lifts the energy degeneracy between states of same spin 
associated to ±t and affects the matching conditions so 
that asymmetry arises. These effects can be seen as chi¬ 
rality phenomena for scattering, analog to magnetic cir¬ 
cular dichroism for the optical absorption in ferromag- 
nets. It makes understandable why the association of 
propagating and evanescent wavevector components may 
exhalt spinorbitronic effects. 

We have presented theoretical evidence for large in¬ 
terfacial scattering asymmetry of carriers vs. incidence 
in semiconducting exchange steps and tunnel barriers. 
The effect appears to be robust. Direct experimental in¬ 
vestigations can probably be performed through angle- 
resolved photoemission spectroscopy. After averaging 
over incoming states, a current perpendicular to the bar¬ 
rier is significantly deflected upon tunneling resulting in 
Tunnel Hall effects and paving the way to original func¬ 
tionalities. Preliminary results indicate that the tunnel¬ 
ing asymmetry in the valence band may play an impor¬ 
tant role in the analysis of tunneling anisotropy magne¬ 
toresistance data m- 
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